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Abstract Neutral ytterbium (Ybl) and singly ionized 
ytterbium (Ybll) is widely used in experiments in quan- 
tum optics, metrology and quantum information sci- 
ence. We report on the investigation of isotope selective 
two-photon-ionization of Ybl that allows for efficient 
loading of ion traps with Ybll. Results are presented 
on two-colour (399 nm and 369 nm) and single-colour 
(399 nm) photoionization and their efficiency is com- 
pared to electron impact ionization. Nearly determin- 
istic loading of a desired number of Ybll ions into a 
linear Paul trap is demonstrated. 



1 Introduction 

Atomic and singly ionized ytterbium is the element of 
choice in numerous experiments in atomic physics and 
quantum optics. Trapping and laser cooling of 1 Yb + 
in Paul traps was initially motivated by its potential 
use as a frequency standard in the microwave or opti- 
cal regime - In addition, this element 
is now employed by an increasing number of research 
groups for implementing experiments related to quan- 
tum information science (e.g., HUJHI1II11IIMI1III1JII3 
[T7]). Neutral Yb is used, for instance, in Bose-Einstein- 
condensation [T5J , photoassociation [19] . or for experi- 
ments aiming at the measurement of parity violation in 
atoms [20] . 
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An important step in trapping ions is the ioniza- 
tion process. For almost 25 years since the first suc- 
cessful trapping and laser cooling of atomic ions [2T1122] 
electron impact ionization was the process of choice. 
In recent years photoionization has begun to estab- 
lish itself as an attractive ionization process for differ- 
ent elements, as it offers distinct advantages over elec- 
tron impact ionization [2511211251126112711281113] . First, it 
is highly efficient in comparison with electron impact 
ionization, typically yielding loading rates that are sev- 
eral orders of magnitude higher. This allows to reduce 
the neutral atom flux considerably, thus avoiding con- 
tamination by deposition of atoms on trap electrodes or 
other, isolating trap components. Second, the electron 
bombardment introduces excess charge into the trap 
volume and on trap components that disturbs the trap- 
ping potential. This effect, too, is avoided when using 
photoionization. And third, by employing a resonantly 
enhanced process, photoionization can be made isotope 
selective in contrast to indiscriminate trap loading with 
the electron impact method. Thus, the need for isotope- 
enriched samples, often a prerequisite for experiments 
with isotopically pure ionic crystals, is eliminated. 

In this article we present the investigation of pho- 
toionization of atomic ytterbium, with particular atten- 
tion to its use as a loading process for ion traps. A fur- 
ther advantage of the scheme presented here is the need 
for only one light source at a convenient wavelength 
near 399 nm that is reasonably close to the now estab- 
lished industrial standard at 405 nm for laser diodes 
making access to inexpensive laser light sources more 
likely. Recent progress in laser diode production pro- 
cesses, however, makes outliers, which deliver the de- 
sired wavelength, less probable and makes it necessary 
to use cooled laser diodes [T4"] . 
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In section [2] we report on laser spectroscopy of the 
x So - 1 Pi resonance near 399 nm in Ybl with the hy- 
perfine structure resolved. Resonance fluorescence spec- 
tra are recorded employing a collimated thermal atomic 
beam emerging from an Yb-oven. We present a simpli- 
fied model of this Yb-oven in section [3] that predicts 
measured oven temperatures satisfyingly and agrees 
with density and flux measurements presented in sec- 
tion |U Resonance enhanced photoionization of Ybl us- 
ing an additional light field near 369 nm is then reported 
in section [5] while in section [S] it is shown that the light 
field near 399 nm alone may ionize Ybl in the presence 
of the rf field used for ion trapping. Section [7] concludes 
this article. 



2 Laser spectroscopy of neutral Yb 

Photoionization of Yb by one photon would require ra- 
diation with wavenumber larger than 50443 cm -1 [29"] 
corresponding to a wavelength below 200 nm in the 
vacuum-uv range. Generating and guiding radiation in 
this wavelength regime is rather inconvenient. In addi- 
tion, one-photon-ionization does not take advantage of 
bound energy levels of Yb in order to make this process 
isotope selective. Therefore it is useful to first excite 
Yb to a bound excited state and subsequently ionize 
the excited atom. 
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Fig. 1 Partial energy level scheme of Ybl. Two resonances con- 
venient for laser excitation and subsequent photoionization are 
indicated by arrows. 



Fig.Q]depicts relevant low lying energy levels in Ybl. 
Starting from the ground state several energy lev- 
els may serve as a first excitation level from which the 
actual ionization process occurs. The inter-combination 
line near 555.6 nm has been extensively studied [301I3T] 
as starting point for a photoionization process. Unfortu- 
nately, this wavelength is presently not accessible using 
an inexpensive diode laser. Therefore, we concentrate 
on the optical dipolc transition from the ground state 
to the : Pi level near 398.9 nm. The x Pi state has a 




lifetime of 5.464 ns [TTJ] and the transition 1 So - x Pi is 
saturated by a light intensity of 60 mW/cm 2 [32] ■ Pos- 
sible decay channels of the x Pi state arc not indicated 
in Fig. [1] as the ratio between the rate for radiative de- 
cay into the ground state and the sum of the rates for 
decay into all other states is larger than 10 7 [53] . 
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Fig. 2 Two-photon- ionization of Ybl via the Pi state. 

As indicated in Fig. [5] a second light-field with wave- 
length below 394 nm is required to reach the continuum. 
This process has been proposed by Sankari et al. [34] . 
Fortunately, this second light-field is easily supplied by 
the cooling laser for the Ybll ion near 369 nm which is 
already present at experiments working with trapped 
Ybll. This makes this scheme very practical indeed, 
as only one other diode laser needs to be set up. This 
process is termed two-colour-ionization and its experi- 
mental investigation is described in detail in section [S] 

This contrasts with the scheme presented in sec- 
tion [TJ] that we label one-colour- ionization. As is shown 
experimentally, a second light- field near 369 nm is not 
necessary to ionize Ybl atoms in an ion trap, instead 
employing only the laser near 398.9 nm is sufficient. 

In the remainder of this section we will be concerned 
with the spectroscopy and isotope selective excitation 
of the 1 So - 1 Pi transition in Ybl. 

Laser induced resonance fluorescence spectra of Ybl 
are recorded using a thermal atomic beam emitted from 
three different ovens heated resistively. The ovens are 
molybdenum tubes with an inner diameter of 0.78 mm 
and a wall thickness of 0.2 mm, spot- welded to a tan- 
talum wire. The ground connection is cut from a tan- 
talum foil (thickness 0.2 mm) with a size of approxi- 
mately 2x4 mm 2 , spot-welded close to the opening of 
the oven tube. Two ovens contain samples of Ybl that 
are isotopically enriched with 171 Yb and 172 Yb , re- 
spectively while the third oven contains a sample with 
the natural abundance of isotopes of ytterbium. Two 
different experimental set-ups are used for the experi- 
ments reported in this article: The isotopically enriched 
ovens are mounted in a vacuum chamber maintained at 
a pressure below 10 _ mbar (when the oven is turned 
off) that also houses a miniature Paul trap (compare 
section I5.II) . The third oven is contained in a vacuum 
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chamber at a typical pressure of 5 ■ 10~ 9 mbar together 
with a linear Paul trap [T3] . 

For excitation of the - 1 Pi transition a tuneable 
External Cavity Diode Laser is used with an emission 
bandwidth of about 2ir x 500 kHz. The wavelength of the 
laser radiation is measured using a wave meter to a pre- 
cision of 0.05 pm using a temperature stabilised HeNe 
laser as a reference corresponding to an uncertainty of 
«95MHz in absolute frequency. The laser beam crosses 
the atomic beam at an angle close to 45° or 90°. 

Resonance fluorescence of Ybl close to 399 nm is col- 
lected, imaged, and detected using the setups that are 
in use for photon counting the resonance fluorescence 
close to 369nm originating from the P1/2 - S1/2 transi- 
tion in Ybll. With the setup containing the isotopically 
enriched Ybl ovens light is collected using an objective 
with numerical aperture of approximately 0.4 and then 
imaged onto a photomultiplicr tube (Burlc C31034A2) 
operated in photon-counting mode. The overall detec- 
tion efficiency is about 1%. The setup with the oven 
containing ytterbium with natural abundance of iso- 
topes uses optical elements for light collection with a 
numerical aperture of approximately 0.5 and a photo- 
multiplier (Hamamatsu R5600P) with an overall calcu- 
lated efficiency of 1.3% for the detection of light emitted 
by Ybl and Ybll. 

2.1 Isotope selective excitation of Ybl 

Figure [3J shows laser induced resonance fluorescence 
spectra of the 1 So - 1 Pi transition in Ybl. These were 
measured employing isotopically enriched samples of 
171 Yb or 172 Yb . The abscissa indicates the frequency 
relative to the centre of the resonance corresponding to 
174 Ybl that was found at a vacuum wavelength near 
398.91 nm. 

The excitation line profile in Fig. [3Ji) is Dopplcr- 
broadened and shifted while in Figs. [3J b) and c) 
Dopplcr-broadcning is strongly reduced, with the mea- 
sured lincwidth of r 172 = 2ir ■ 71(1) MHz and A 71 = 
2-7T • 89(1) MHz a factor of approximately 3 higher than 
the natural linewidth of 2ir ■ 28 MHz due to saturation 
broadening and the Doppler-broadening due to the di- 
vergence of the atomic beam. The hyperfme splitting 
of the : Pi level of the 171 Yb isotope is clearly resolved. 
Small peaks to the left and right of the respective main 
resonance in Figs. [3] b) and c) are due to impurities 
from other isotopes. 

After having obtained sufficient spectral resolution 
with the isotopically pure 171 Yb and 172 Yb ovens, a new 
oven with a natural abundance isotope distribution was 
built, and Fig. @] shows the atomic resonance fluores- 
cence spectrum measured with this oven. The spectrum 
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Fig. 3 Laser induced resonance fluorescence spectra of the 1 So - 
x Pl transition in Ybl. An isotopically enriched sample of 171 Ybl 
(a, b) or 172 Ybl (c) was used. The detuning of the laser frequency 
is given relative to the resonance of 174 Yb. a) 171 Ybl: The an- 
gle between laser beam and atomic beam is about 45 degrees. 
The resonance is shifted and broadened by the Doppler effect, 
b) 171 Ybl: The laser beam crosses the atomic beam at right angle 
allowing for resolution of the hyperfme structure. The rest abun- 
dance of 172 YbI can be seen as a small resonance peak to the left 
of the two prominent hyper-fine resonances of 171 Ybl. c) 172 YbI: 
Dopplcr-reduced spectrum. 



was fitted to a sum of six Voigt profiles to obtain the iso- 
tope shift of the resonance frequencies and agreement 
within mutual errorbars was found with |35j . In this 
spectrum the different isotopes are resolved, allowing 
for selective ionization and hence trapping of different 
isotopes of ytterbium. 

The resonance 1 So - 1 Pi , F=3/2 corresponding to 
the isotope 173 YbI is separated from the 1 Sq - 1 Pi res- 



onance in 7 Ybl by only 17 MHz 



Thus it is diffi- 



cult to spectrally resolve these two isotopes. However, if 
loading an ion trap with 172 YbII is desired, then this is 
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Fig. 4 Laser induced resonance fluorescence spectra of the 1 So - 
J Pl transition in Ybl with reduced Doppler-broadcning from an 
oven with natural abundance isotope distribution. All hyperfine 
information displayed at the resonances concerns the upper levels 
of the transitions since the ground state exhibits no hyperfine 
splitting. Nonlinear regression was performed by assuming fixed 
spacing given by 35 and fitting a common horizontal and vertical 
displacement, a natural line width, and residual Doppler width, 
and individual peak heights. 



not an obstacle, since the odd isotope's hyperfine struc- 
ture prevents this ion from being efficiently laser cooled: 
If no microwave radiation is present that avoids opti- 
cal pumping into a 'dark' hyperfine level (with respect 
to laser light near 369nm driving the S1/2 - P1/2 tran- 
sition in Ybll) of the electronic ground state, then it 
is unlikely that 173 Yb + at a thermal velocity is suffi- 
ciently laser cooled to be trapped. Furthermore, these 
ions would be visible as dark gaps in linear ion crystals. 
This has not been observed during the loading of ion 
crystals reported in section [5?3l 



3 One dimensional heat model of the oven 

We have manufactured a variety of ovens with different 
geometries which all have in common that ytterbium 
inside a metal tube with one open end is evaporated by 
ohmic heating. 

The latest ovens were designed to meet the spe- 
cial demand of a micro structured segmented linear 
ion trap. To avoid coverage of the trap chip with neu- 
tral ytterbium, the atomic beam should be restricted 
to a region narrower than the electrode separation. Ini- 
tial experiments with squeezed oven nozzles gave good 
collimation results, confirmed by analysing transverse 
Doppler-broadening profiles. But within 30 hours of op- 
eration the Doppler width broadened gradually, most 
likely due to clogging or re-deposition in the nozzle. 
Satisfying long term behaviour could be achieved by a 



laser cut slit aperture a few millimetres apart from the 
oven nozzle. 



molybdenum wire, = 0.7 mm 

steel wire, locally thinned 




steel tube, 0outer = 1 mm 0inner = 0.6 mm 
Fig. 5 Sketch of the oven 



Our recent ovens are made of a steel tube and wire, 
which is thinned at one point to give a high local ohmic 
loss. A steel tube, partially filled with a small amount 
of ytterbium is mounted close to this point and a col- 
limatcd ytterbium beam emerges from the open tube 
end. The feed wires are made of molybdenum. 

A simple heat model for a small piece of wire states 
that the net flux of power per length element vanishes 
for a steady state temperature distribution 



dPohm - dP Ia d - dP cond = 



(1) 



with the infinitesimal power contributions dP due to 
ohmic heating, radiation, and conduction. These powers 
for a small length interval dx of a cylindrical wire are 
given as 



dPohm = PdR = 



dx 



dP rRd = dAaeT 4 = 2irrae (T 4 - T n 4 ) dx 



dP c 



cond 



: dP cond {x - dx/2) - dP cond (x - dx/2) 

- Kirr 2 (T'(x + dx/2) - T'(x - dx/2)) 

- nirr 2 T" (x)dx 



(2) 
(3) 



(4) 



In the cq. ©-(H]) the heating current is denoted as 
/, p(x) is the electrical resistivity, r is the radius of the 
wire, a the Stcfan-Boltzmann-constant, e the emissiv- 
ity of the wire material and Tq the ambient tempera- 
ture, and the temperature conductivity is termed as k. 
Inserting eq. (J2H|) into ([T]) one finds that the second 
derivative of the temperature has its origin in ohmic 
heating and radiative losses: 



— KTrr 2 T" (x) + 



I 2 p(x) 



- 2nr<7€ (T 4 - T n 4 ) = 



(5) 



When the oven is modelled as a one dimensional system 
made of the feed through wire with a small incision as 
the heating element the electric resistivity along the 
wire is given by 



p(x) = po + Rirr 2 5(x) 



(G) 
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where po is the resistivity of the wire, and R the resis- 
tance of the incision at x = (treated as point like). 
We integrate cq. ([5]) numerically, which allows us to ac- 
count for the temperature dependence of electric and 
heat conductivity. We take into account two boundary 
conditions: i) the temperature at holders at x = ±£ 
equals the ambient temperature To, and ii) the temper- 
ature distribution has a kink at the heating element at 
x = and the derivative changes by I 2 R/Kirr 2 . 
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Fig. 6 Temperature distribution along a wire with an incision 
at x = and held at room temperature at x = ±10 cm 



Figure [6] shows a numerical solution for the temper- 
ature distribution using a current of 1.5 A and parame- 
ters that describe our oven adequately. These parame- 
ters are either obtained from literature (as for k and p) 
or from our own measurements (as for R) . In addition, 
we take into account the linear temperature coefficient 
of the electric resistivity. The parameters used for cal- 
culating and plotting the temperature distribution are: 
r = 0.35 mm, I = 10 cm, R = 0.06 ft, n = 138W/mK 
[36], p = 5.57 -10- 8 (1 + 4.82 -10~ 3 {T -T )/K)mn [37]. 
Figure [6] shows the temperature distribution for various 
linear temperature dependencies of the heat conductiv- 
ity. The distribution is found to be dominantly triangu- 
lar, with negligible ohmic losses in the wire and small 
positive curvature for higher temperatures regions due 
to radiative losses. The model gives reasonable agree- 
ment with measurements of the oven tip temperature 
considering that we did not take into account the ge- 
ometric details of our oven. Best agreement is reached 
when a temperature dependence of the heat conductiv- 
ity is assumed that is approximately a factor 6 lower 
than the coefficient for the resistivity (as it is for iron). 

The oven tip temperature was measured close to the 
incision shown in Fig. fusing a j-type thermistor as a 
function of heating current (see Fig. [7]) and a hyperbolic 



dependence is found. This functional dependence is un- 
expected at first sight: one would expect a quadratic 
behaviour for small currents, as the heating powers de- 
pends quadratically on the current, and radiative losses 
are negligible. In fact, such a dependence is found for 
small currents. For high currents one would expect a 
square root dependence of the tip temperature, arguing 
that dominating radiative losses which scale as T 4 are 
balanced by a heating power behaving as I 2 but such a 
regime is not observed for the currents that we applied. 
An extended linear regime can be explained, however, 
taking into account the temperature distribution along 
the wire as above. The hyperbolic dependence is repro- 
duced by a simple analytical model assuming a trian- 
gular temperature distribution in accordance with the 
results shown in Fig. [51 



T(x) = T n 



T, 



To, 



\x\ < 



(7) 



with the ambient temperature To. The heat loss 
through the wire-ends at x = £ equals 2KTrr 2 T' and 
power balance requires 

I 2 R = 2nnr 2 T' + [ 27rrcre (T(a;) 4 - T 4 ) dx (8) 
J -l 

(9) 



= 2Kirr 2 T' 



Aitraei 



5(T n 



To) 



(^max ^0 ) 



After polynome division of the radiative term and con- 
sidering terms up to second order in T max the peak 
temperature is obtained as 




T n 



with the constants a and /3 given by 
R 



a 



/3 



(10) 



(11) 



which gives the observed hyperbolic dependence. We 
have measured the temperature of the oven tip with 
a vacuum compatible thermistor and find good agree- 
ment with a hyperbolic model. For a current above 1 A 
(that leads to a useful atom flux as will be reported 
below), a linear relation between current and temper- 
ature is found. It should be noted that the measured 
quantitative asymptotic behaviour is not predicted by 
the coefficients extracted from (fTTj). most likely due to 
geometric details and material properties that are ne- 
glected or simplified in our model. However, the im- 
portant general insight remains that heat conductivity 
and radiative losses yield a polynomial expression in 
T max — T) in cq. ([5]) which results in a hyperbolic de- 
pendence of the maximum temperature as a function of 
heating current. 
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Fig. 7 Temperature of the oven tip as a function of heating 
current. 



4 Neutral Atom Densities and Atomic Flux 

The temperature of ytterbium in a closed oven would 
result in a vapor pressure, which, for a temperature 
independent latent heat, is predicted by the Clausius- 
Clapeyron equation to take the simple form 



Pvp(T) = Po e- Tl/T « 353334 Pa- e ^ 9676K / T 



(12) 



where the numerical constants on the right-hand-side 
were obtained by to match measured vapor pressures 
in the temperature range 623 K to 931 K [35]. We used 
this relation despite the extended temperature range 
investigated here when modelling atomic density and 
flux data. Assuming the ytterbium vapor to behave as 
an ideal gas, the particle density can be calculated from 
the vapor pressure 



Pv P (T) _p e 



-Tx/T 



(13) 



Furthermore, for Maxwell distributed atom velocities, 
one can calculate the flux of atoms through an circular 
aperture with radius r 



8rrr 4 
mkvtT 



Po e 



-Tx/T 



knT 



[v) po e 



-Ti/T 



with the mean velocity of the atoms in the beam 



(v) = 



8k B T 



(14) 



(15) 



When calculating densities and fluxes from currents us- 
ing equations (fTUf and (|12l - [Ft]) , one has to bear in 
mind, that the effective temperature might in general 
be lower than the one obtained for the heating element, 
since the atoms can thermalize with the whole, and gen- 
erally cooler, oven tube. 



Atom densities and fluxes were characterized for 
the oven which was optimized for use with a micro- 
structured segmented linear Paul trap. The neutral 
atom fluorescence is imaged onto a small aperture, 
which in turn is imaged onto a photomultiplier. Mag- 
nification and aperture are chosen such that the inten- 
sity and atom density can be assumed approximately 
constant over the detected area. We take into account 
saturation by determining experimentally the intensity 
profile of the laser and the total power. For all measure- 
ments we assume a steady state temperature distribu- 
tion which is ensured by a 10 minute warm-up. 

The measured current dependence (Fig. [5J and 
of the density and flux is fitted by a simple model: 
the atomic vapour is in thermal equilibrium with the 
oven. The oven has an effective temperature, which 
is assumed to exhibit a linear dependence on current 
as found with thermistor measurements for higher cur- 
rents (compare Fig. [7]). The dependence of the effective 
oven temperature is not necessarily identical with the 
dependence of the oven tip temperature: we allowed 
for a temperature offset and a different slope; in gen- 
eral, the effective temperature is expected to be lower 
than the one measured by the thermistor. The tem- 
perature is taken to calculate the vapour pressure and 
flux. Both quantities are multiplied by a factor, taking 
into account the finite atomic beam divergence, non- 
equilibrium effects, imperfect characterization of our 
detection efficiency and saturation. 

The model describes the measurements quite well. 
It states that the effective temperature of the oven (and 
thus the atoms' temperature) shows the same slope 
with respect to the current, but is substantially lower 
than the temperature of the oven tip (by approximately 
50K). 



5 Two-color-ionization of Ybl 

In this section we report on resonance enhanced isotope 
selective ionization of Ybl with subsequent loading of 
the ions into a Paul trap (ring trap or linear trap). The 
possibility to selectively excite a certain isotope out of a 
mixture of natural abundance material is an advantage 
of photoionization. This selectivity is achieved by the 
first excitation step from the ground state of Ybl into 
the 1 Pi state as reported in the previous section. The 
second excitation step from the : Pi state of Ybl into the 
ionization continuum is achieved by means of laser light 
near 369 nm (Fig. [5]). The dependence of trap loading 
efficiency on laser intensities, laser frequencies, and the 
relative flux of neutral atoms is studied. 



7 




1.8 1.9 2 2.1 2.2 

oven current (A) 



Fig. 8 The atom density from a collimated ytterbium with nat- 
ural isotope abundance as a function of oven heating current. 
The beam is collimated by an aperture separated 3 mm from the 
oven nozzle. The density was measured approximately 5 mm away 
from the collimating aperture. 
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oven current (A) 



Fig. 9 The atomic flux density from a collimated ytterbium with 
natural isotope abundance as a function of oven heating current. 
The density data presented in Fig. [5] was combined with measure- 
ments of the line shift due to the mean beam velocity as shown 
in Fig.[3ji). 

5.1 Experimental procedure 

A part of the experimental setup is described already 
in the preceding section that deals with resonant ex- 
citation of Ybl. Now, taking advantage of this isotope 
selective excitation of the x Pi state, we are concerned 
with photoionization of Ybl and trapping of Ybll. Ions 
are trapped in either a miniature ring Paul trap or a lin- 
ear Paul trap. The ring trap electrodes are made from 
molybdenum wire with a diameter of 0.3 mm. The ring 
electrode has an inner diameter of 2.0 mm and the end- 
cap electrodes are 1.4 mm apart. An rf field at 9.5 MHz 
with peak amplitude of about 700 V is applied between 
ring and endcap electrodes. The linear trap is made of 



four Molybdenum rod electrodes 0.5 mm in diameter, 
held by two Macor plates spaced 6 mm apart. The ra- 
dial distance from the axis to the surface of the rods is 
0.75 mm. Two Molybdenum end-caps with a diameter 
of 0.4 mm serve for axial confinement, and are spaced 
4.1mm apart. An rf field with frequency between 10 
and 20 MHz and amplitude of up to 2 kV is applied to 
the rod electrodes for radial confinement, and a low 
voltage of several volts only is applied for axial confine- 
ment. Typical secular frequencies are 30-60 kHz in the 
axial direction and 600-800 kHz radially. 

Light near 369nm driving the S1/2 - P1/2 transi- 
tion in Ybll is supplied by a frequency doubled com- 
mercial Ti:Sa laser (Coherent MBR110). This will be 
termed "the cooling laser" in the remainder of this ar- 
ticle. A diode laser delivers light near 935 nm and drives 
the D3/2 - [3/2] 1/2 transition in Ybll to avoid optical 
pumping into the metastablc D 3 / 2 state [3"9"] . 

An electric current resistively heats the oven, and 
after some time, a stationary temperature is obtained. 
Therefore each oven was heated for at least 10 min prior 
to the actual experiment, unless stated otherwise. The 
oven temperature determines the vapour pressure and 
atom density inside the oven and the atomic flux den- 
sity in the particle beam that emerges from the open 
end of the oven. 

In order to determine quantitatively the atomic 
number density in the trapping region of the ion trap 
one would need knowledge of i) the spatial profile of the 
atomic beam and the laser beam in the region where 
they cross, ii) the overall detection efficiency of the sys- 
tem measuring the resonance fluorescence in this cross- 
ing region as function of the position of the fluorescing 
atoms (the depth of focus for the optics collecting the 
fluorescence is rather small due to a large numerical 
aperture), and iii) the transition line strength of the 
1 So - x Pi transition. Thus, deducing the number den- 
sity from the observed resonance fluorescence spectra 
is in principle feasible, albeit with a large systematic 
error. If the number density were known, the investi- 
gation of ion loading rates as a function of oven cur- 
rent reported below could be presented as a function 
of neutral atom density instead. The ion loading rates 
depend, in addition on i) the ionization cross section 
for the different ionization schemes reported (one-color, 
two-color, electron impact), ii) the quasi static electric 
field experienced by the atoms in the trap region that 
varies spatially and in time in a different manner for 
each trap, and iii) the rate at which different Ybll ions 
are laser cooled which in turn is a function of the spatial 
profiles and overlap of the laser beams near 369 nm and 
935 nm used for laser cooling Ybll and their frequency, 
intensity, and polarization (in addition, for odd isotopes 
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also of the parameters characterizing the microwave ra- 
diation necessary to avoid optical pumping) . We do not 
attempt here to give a full quantitative description that 
might be plagued by a large uncertainty while not al- 
tering the conclusions drawn from the measurements 
reported below. 

The wavelength of the laser near 398.9 ran was set to 
a value corresponding to the 1 So - 1 Pi resonance of the 
desired isotope to be loaded into the trap, and the peak 
laser intensity (Gaussian beam profile) in the ionization 
region was 3 W/cm 2 , if not explicitly stated otherwise. 
For brevity we will label this laser "excitation laser" in 
what follows. 

When both the laser near 369 nm and near 398.9 nm 
are directed into the trap region the atomic resonance 
fluorescence signal is added on top of the ionic fluo- 
rescence. Experimentally the signals due to Ybl and 
Ybll are discernible by blocking temporarily the re- 
pump laser at 935 nm which leads to optical pumping 
into the metastable D 3 / 2 state of Ybll and thus the 
ionic fluorescence is interrupted. Typically, the fluores- 
cence signal collected from neutral atoms amounts to - 
depending on oven current - several tens of kHz. 



5.2 Two-color ionization 




10 20 30 40 50 

time (s) 



Fig. 10 Two-color photoionization of 17 Ybl as a function of 
time for different electric currents heating the atom oven. The 
fluorescence signal is proportional to the number of trapped Ybll 
ions for a low number of ions, i.e., short loading times (here, a 
single ion produces a photon count rate of xs 15 kHz). In this 
region the signals are fitted by a linear function to extract loading 
rates as a function of oven current. Steps in the photon count 
rate for small heating currents indicate an increase by one of 
the number of trapped ions. The oven currents were (in order of 
ascending signal slope): 1.8 A, 1.9 A, 2.0 A, 2.18A, 2.3 A, 2.47 A, 
2.57 A. 



Fig. [TU] shows combined atomic (Ybl) and ionic 
(Ybll) fluorescence signals in dependence on the ioniza- 
tion time and the oven current for two-color-ionization. 
The procedure to record the time series in Fig. [10] is as 
follows: The ion trap is emptied by turning off the rf 
trapping field and all lasers except the excitation laser 
at 398.9 nm are turned on. After 3 s the excitation laser 
is unblocked, leading to fast loading of the trap and 
accordingly to an increase in the ionic fluorescence sig- 
nal. With a fluorescence rate of 15 kHz for a single ion 
here, on the order of a 100 ions are trapped in the first 
15 seconds for the highest oven current at 2.57 A, in- 
dicating the high efficiency of this process. This oven 
current represents the usual value used with electron 
impact ionization. 

The resonance fluorescence signal begins to satu- 
rate, after a time depending on oven current, generally 
for ion numbers greater than about 20, owing both to 
the limited trapping capacity and the reduced efficiency 
of the optics for imaging a large ion " cloud" (the optical 
set up used here is optimized for imaging a single ion 
located in the center of the trap). A limited capacity 
for trapping can be attributed to the laser beam waist 
being too small to illuminate large clouds. In addition, 
large hot clouds exhibit a large number of micromo- 
tion sidebands, such that the scattering of resonance 
fluorescence induced by a laser with narrow emission 
bandwidth is reduced dramatically. 

The oven current is varied between 1.80 A and 
2.57 A. For those currents we obtained feasible loading 
rates and low vacuum contamination. As expected, the 
rate of loading ions increases with the atom flux den- 
sity. Depending on the probability for ionizing a neutral 
atom passing by, two dependencies are conceivable: for 
almost unit probability of ionization, the loading rate 
will be proportional to the atom flux, whereas for low 
probability of ionization, the loading rate can be ex- 
pected to behave proportional to the atom density. A 
comparison of loading rates on the order of a few per 
second and much higher measured fluxes shows that the 
latter is the appropriate model, and Fig. [12] shows the 
measured loading rate together with a fit that assumes 
proportionality to the atom density. 

In Fig. [TT] the loading rate using two-color- 
photoionization is compared with electron impact ion- 
ization. The time series indicating the fluorescence sig- 
nal when using electron impact ionization shows a sharp 
increase at about 142 s that signifies the trapping of 
a single 172 Yb ion. In contrast, the ionic signal rises 
quickly when using photoionization indicating trap- 
ping of more than 100 ions during the first seconds 
after turning on the excitation laser. The two-color- 
ionization process results in an estimated increase in 
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50 100 150 200 

time (s) 

Fig. 11 Comparison of efficiency for electron impact and two- 
color-ionization of 172 YbI (oven current set to 2.5 A) 1131 . The 
upper trace (dots) shows the rise in fluorescence when two color 
photoionization is used and a fit (solid line) to the linear slope of 
the signal as in Fig. 1 101 (note that the linear slope is not displayed 
as a straight line due to the logarithmic depiction). The step 
in the lower time trace (electron impact ionization) signifies the 
trapping of a single ion. 

trap loading rate of about 3000 in comparison with elec- 
tron impact ionization. For this measurement the oven 
current was set to 2.5 A. 

effective temperature (K) 




oven current (A) 

Fig. 12 Loading rate extracted from the slope of fluorescence 
increase in the linear regime as a function of oven current (com- 
pare Fig. 11011 . A thermodynamic vapor density is fitted to the 
data assuming a linear dependence between current and effec- 
tive temperature as observed from probe measurements (see sec- 
tion^. The fitted dependency is comparable to the one measured 
in section[4] even though the oven was constructed with a slightly 
different design. 

In the two paragraphs to follow we examine the in- 
fluence of the frequency of both the laser light field 
near 369 nm and the one near 399 nm on the efficiency 
of the ionization process. The ionization model depicted 
in Fig. [2] is corroborated by these measurements. 



Laser Detuning A (MHz) 




£ -600 -300 300 600 £ 

Laser Detuning A (MHz) 



Fig. 13 Dependence of ionic resonance fluorescence signal - and 
hence the loading rate - on detuning of both lasers for two-colour- 
photoionization during the loading time. The fluorescence signal 
shows a strong dependence on the detuning of the 398 nm laser 
(filled dots), since this excitation is a resonant process. The solid 
line indicates a fit using a Voigt profile. The non-resonant exci- 
tation with the 369 nm laser above ionization threshold shows no 
dependence upon variation of wavelength (open circles) and can 
be fitted by a constant (dashed line). The ionization time was 
30 s at 2.59 A oven current. 

Fig. [T2] shows the resonance fluorescence signal of 
trapped 172 Yb ions under variation of the frequency 
of the cooling laser close to 369 nm over a range of 
800 MHz during the loading phase. The scan was taken 
with the oven current set to 2.59 A and 30 s ionization 
time. During the observation stage of the trapped ions, 
the laser at 369 nm was set close to resonance. The ob- 
served fluorescence signal is consistent with a structure- 
less flat dependence of the ionization rate on the detun- 
ing of the cooling laser in the frequency range depicted 
here. This is in accord with the ionization model that 
the atom is excited from the 1 Pi state into the ion- 
ization continuum (Fig. [2]). The energy of the photon 
at 369 nm is larger than the energy needed to reach 
the continuum (equivalent to a wavelength of 394 nm) , 
thus no dependence on small frequency variations is ex- 
pected. 

On the other hand, based on the excitation spec- 
tra of neutral Ybl the detuning of the excitation laser 
at 398.9 nm is expected to have a distinct influence on 
the ionization efficiency for a given isotope. The second 
data set in Fig. [T3] confirms this assumption. There, 
the excitation laser is scanned over a range of 1.2 GHz 
and the resulting ionic fluorescence is recorded. This 
graph may also be compared to Fig. [3k,, since both data 
sets were taken using the 171 Yb isotope and an angle 
of 45 deg between the atomic and laser beam, respec- 
tively (thus the line profile is dominated by Dopplcr- 
broadening). The oven current was set to 2.82 A and 
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the time during which the excitation laser was turned 
on was 60s. 

5.3 Nearly deterministic loading of an ion trap 

The last paragraph of this section deals with an exper- 
imentally important feature of photoionization of Ybl 
for ion trap experiments: It is desirable to (nearly) de- 
tcrministically load a desired number of ions into the 
trap. In addition, the flux of neutral atoms possibly 
contaminating trap structures should be significantly 
reduced. So far, the oven current was set to a certain 
value and kept there to allow for comparison of dif- 
ferent loading rates. Now, when trying to quickly load 
ions, all laser frequencies and intensities are set first 
before the oven current is turned on. The wavelength 
of the excitation laser near 399 nm is set with an ac- 
curacy better than 80 MHz using a wave meter, which 
is sufficient to find the atomic resonance without hav- 
ing to rely on the atomic beam. Thus, already the first 
atoms emitted from the oven have a chance of being 
ionized by efficient two-color-ionization and are imme- 
diately laser cooled by the same laser light that ionizes 
them. It turns out that the loading rate using this pro- 
cedure is such that any desired (low) number of ions is 
easily selected simply by blocking the ionization laser 
at the right time. 




o ' 1 1 1 ' 

50 100 150 200 

t(s) 

Fig. 14 Deterministic loading of a desired number of ions by 
two-color-photoionization 1131 . The loading rate increases with 
time, since the oven emitting Ybl atoms heats up during the 
measurement 

Fig. [14] shows the number of trapped ions plotted 
against time. This experiment was performed with a lin- 
ear Paul trap and the ionic fluorescence was detected by 
an image-intensified CCD-camera. The time before the 
first ion is detected varies between 60 - 180 s after the 



oven-current is turned on, depending on laser parame- 
ters. Then, while the oven is being heated the loading 
rate increases. This graph shows that nearly determin- 
istic loading of a desired number of ions is possible using 
this procedure. 

6 Single-color-ionization of Ybl 

The energy carried by two photons at 398.9 nm is not 
sufficient to bridge the energy gap between the ground 
state of Ybl and the ionization continuum. Therefore, 
it came initially as a surprise that Ybll ions could be 
produced using this laser light alone. 

For the experiments presented in this section all 
lasers are blocked first for a prescribed time while 
the laser at 399 nm is on. Then the laser exciting the 
atomic resonance is blocked while lasers near 935 nm 
and 369 nm are turned on for cooling and detecting the 
ion. 

Comparing the two methods shows that for the 
experimental parameters used here, the two-color- 
ionization is almost two orders of magnitude more ef- 
fective than the one-color-ionization. In the latter case, 
though, the ions are cooled only after the ionization pro- 
cess is finished when the cooling lasers are unblocked. 
For the two-color scheme, on the other hand, the ions 
are cooled as soon as they are ionized. 

As an explanation of why ions are produced by 
399 nm light alone, 3-photon-ionization would be con- 
ceivable, with the first photon exciting the atomic 
- 1 Pi resonance as before, and two more photons ioniz- 
ing the atom. This "second" two-photon process could 
again be resonance enhanced, if the 399 nm light excites 
a Rydberg state in Ybl. If the absorption of the last 
two photons is not resonance enhanced, then a rough 
estimate shows that at a light intensity of a few Watts 
per cm 2 this three-photon-ionization rate (proceeding 
resonantly with the first photon through the : Pi state 
and then with two more photons non resonantly into 
the continuum) should be about 12 orders of magni- 
tude lower than the two- photon rate (i.e., proceeding 
resonantly with the first photon through the x Pi state 
and then with one more photon into the continuum) 
|40j . Since such a low ionization rate is not consistent 
with our experimental observation, we conclude that a 
three-photon process that is not a second time reso- 
nance enhanced cannot be responsible for the observed 
production of Ybll using light near 399 nm. 

Is there a Rydberg state that could be resonantly 
excited starting from the 1 Pi state with the same 
light that excites the 1 So - 1 Pi transition, and thus 
lead to a three photon process that is resonantly en- 
hanced a second time? From this Rydberg state, a 
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third photon at the same wavelength would then ion- 
ize the atom. Starting from the 1 Pi level the dipolc- 
selection rules allow transitions into either S- or D- 
Rydberg states. Checking the energy levels of the S- 
and D-Rydberg series shows that no level exists at 
50136.444 cm -1 , the energy corresponding to two pho- 
tons near 398.9 nm j4Hl4"2"] . The nearest unshifted sin- 
glet states are the 6s23s 1 So at 50130.98 cm -1 and the 
6s22d 1 D2 at 50148.59 cm -1 , with an energy difference 
to the value above of 5.46 cm -1 and 12.15 cm -1 , respec- 
tively. We therefore exclude resonant excitation of an 
unshifted Rydberg state as a likely candidate for the 
observed photoionization by a one-color light-field. 

When considering possible ionization pathways we 
have neglected so far the influence of the rf trapping 
field on the ionization process. Given the trap dimen- 
sion of the ring trap and the applied rf voltage the peak 
electric field is estimated as 7 x 10 5 V/m. With the 
Stark shift taking on a value of about 10 -6 Hz/ (V/m) 2 
for Hydrogen with n = 2 and being proportional to 
n 7 (with principal quantum number n) |43j it is ex- 
pected that the Stark shift of the 1 Pi state is negligible 
here while Rydberg states with n ~ 22 arc considerably 
shifted in energy. In fact, the trapping potential lowers 
the ionization threshold such that a second photon at 
398.9 nm directly reaches the continuum. 

In a simplified picture a Hydrogen atom exposed to 
a constant external electric field in direction x shows an 
asymmetric binding potential along this direction as 

e 2 1 

V = -— eEx, (16) 

47re x K ' 

with elementary charge e, electron nucleus separation x, 
and the magnitude of the external electric field E. The 
potential maximum in the direction of the electric field 
is lowered with respect to the unperturbed atom which 
corresponds to a lowering of the ionization threshold by 



with the peak voltage Uq and the electrode separation 
2ro yielding a peak electric field of 



AU = 



7re ' 



(17) 



This is reminiscent of the Schottky effect [44] which 
describes the lowering of the work function of solid state 
materials in the presence of an electric field. 

The difference of the ionization energy U = 
50443 cm -1 and the energy of two photons 2E 39S = 
50136. 4cm~ 1 requires a lowering of the ionization en- 
ergy by at least AU = 306.6cm -1 . Using eq. (|T7)) . 
this corresponds to a required field strength exceeding 
|E ph | = 2.5- 10 5 V/m. 

On the other hand the peak potential in a linear ion 
trap along the z direction is given by 



E 



peak 



2 W* 2 + y 2 = 2C/ - (19) 



Thus the radius, at which the required electric 
field strength for one-color-photoionization is exceeded 
amounts to 



|Eph|?- 2 
2C/n ! 



(20) 



^pcak 



(18) 



which gives for our trap (?7o = 1000 V and 
2 tq = 1.5 mm) a separation from the trap center of 
r = 140 pm. This gives .75 eV of initial quasi potential 
energy which has to be laser cooled. 

A full quantitative treatment of the electric field 
enhanced photoionization has to account, in particular, 
for the density and velocity distributions of the Ybl 
atoms moving through the oscillatory trapping field, 
and the exact spatial dependence of the trapping field. 

The significantly lower trap loading rate when us- 
ing 399 nm light alone for photoionization, as compared 
to the two-color process, is explained when considering 
i) that only a fraction of the atoms interacting with 
399 nm light has their ionization threshold sufficiently 
lowered to be ionized, and ii) the fact that during the 
short time interval between turning off the laser near 
399 nm and turning on the lasers necessary for cooling 
Ybll a fraction of the ions will have already left the 
trapping region. 

Fig. [T5] shows the dependence of the loading pro- 
cess as a function of the detuning of the light near 
399 nm from the 1 So - 1 Pi resonance in 172 YbI. One 
dataset (filled circles) was taken with an oven current 
of 2.59 A and the ionization time of 30 s. Due to an 
angle of 45 degrees between laser and atomic beam, 
the line width (FWHM 320 MHz) shows substantial 
Doppler-broadcning (when comparing this with Fig. [3^,) 
it should be noted that (Fig. [3^) shows a resonance flu- 
orescence spectrum of 171 Yb, while the data in Fig. [TS] 
shows the ionization of 172 Yb as a function of laser de- 
tuning) . The second scan was taken with the laser beam 
perpendicular to the atom beam, an oven current of 
2.79 A and ionization time of 30 s. The line width dis- 
played here (FWHM 110 MHz) is consistent with the 
width of the atomic resonance line recorded as shown in 
Fig-Et) with 172 YbI. Thus, the ionization rate as a func- 
tion of frequency of the laser light field near 399 nm is 
consistent with the absorption of the first photon ( 1 So - 
*Pi transition) being a resonant process and the second 
photon absorption ending in the ionization continuum. 
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Fig. 15 Ionization efficiency of one-colour photoionization of 
172 Yb as a function of detuning of the 399 nm light from the 1 So 
- J Pi resonance. The filled circles show the fluorescence signal as 
a function of laser detuning, similar to Fig. 1131 The ionization 
time was 30 s at an oven current of 2.59 A and laser and atomic 
beam enclose an angle of 45 degrees. Fitting a Voigt profile yields 
a FWHM of 320 MHz. 

The open circles display the number of trapped ions after a load- 
ing time interval of 30 s at 2.79 A oven current Here laser and 
atomic beam are perpendicular yielding low Dopplcrbroadcning. 
Fitting a Voigt profile yields a FWHM of 110 MHz which is con- 
sistent with fluorescence profiles from neutral 172 Yb . 

7 Conclusion 

We report on two photoionization processes for neu- 
tral ytterbium that allow for loading of Yb + -ions into 
a Paul-trap. Both methods are implemented using just 
one additional diode laser near 399 nm. Laser diodes 
exist for this wavelength, and are close to the DVD- 
standard at 405 nm. Since the saturation intensity of 
the 1 So - 1 Pi transition in neutral Yb is rather small, 
a diode delivering of the order 10 mW is sufficient for 
efficient ionization. Thus, a usual diode laser system in 
external cavity setup characterized by a small emission 
linewidth and easy to tunc suffices to implement effi- 
cient photoionization of Yb. 

The second step in the photoionization process (af- 
ter excitation of the 1 So - 1 Pi resonance) is achieved by 
absorption of a photon near 369 nm, the same light field 
used for cooling and detection of Ybll ions. Therefore, 
it is ensured that the ionization process takes place in 
a region where the ion is immediately laser cooled and 
thus trapped and detected. 

In addition, it is demonstrated that in the presence 
of the rf field used for trapping ions in a Paul trap, the 
light field near 399 nm alone suffices for a second exci- 
tation step into the ionization continuum shifted by the 
quasi static electric field applied to the trap electrodes. 

The isotope-shift of the 1 So - 1 Pi transition in neu- 
tral Ybl is resolved in excitation spectra when reduc- 



ing Dopplcr-broadcning by crossing the atom and the 
laser beam at an angle near 90°. Thus isotope selec- 
tive ionization is possible. The efficiency of two-color- 
photoionization is found to be larger by a factor of 
3000 in comparison with electron impact ionization for 
the experimental set-up used in these experiments. We 
showed that (nearly) deterministic loading of a desired 
number of Ybll ions into a linear Paul trap is feasible. 
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the manuscript. Financial support by the Deutsche 
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tegrated Project QAP, STREP Microtrap), and by 
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